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ABSTRACT 

The Palomar High-precision Astrometric Search for Exoplanet Systems monitored 51 subarcsecond 
binary systems to evaluate whether tertiary companions as small as Jovian planets orbited either the 
primary or secondary stars, perturbing their otherwise smooth Keplerian motions. Six binaries are 
presented that show evidence of substellar companions orbiting either the primary or secondary star. 
Of these six systems, the likelihoods of two of the detected perturbations to represent real objects are 
considered to be "high confidence" , while the remaining four systems are less certain and will require 
continued observations for confirmation. 

Subject headings: astrometry - binaries: close - binaries:visual - techniques:interferometric 



1. INTRODUCTION 

The use of astrometric measurements to detect the re- 
flex motions of stars caused by substellar companions 
orbiting t hem has a long history filled with false alarms. 
Famously, |van de Kamp ( 1963 1 claimed to have discov- 
ered giant planet companions to Barnard's Star. His 
first estimates of a single planet of 1.6 times the mass of 
Jupiter with an orbital period of 24 years and an eccen- 
tricity of 0.6 were later revised to two planets with orbital 
masses a nd periods of 1.1 Mj at 26 years and 0.8 Mj at 
12 years ( van de Kamp| 1969|. He never accepted grow- 
ing evidence from other astronomer s that these discov- 
eries were not repeatable elsewhere (Gatewood & Eich- 
horn||1973 Hershey|l973 ) ; today it has be en shown coi> 
clusiy ely that these planets are not real (Kiirster et al. 
20031. 



Han et al. (2001) used Hipparcos measurements to an- 
alyze stars with known radial velocity (RV) detected ex- 
oplanet candidates. The precision of Hipparcos was in- 
sufficient to detect the reflex motions if the objects are 
in fact planets; however, if instead the orbits are face-on, 
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the actual companion masses would be larger than the 
RV-derived masses, so the resulting much larger motions 
could have been detected by Hipparcos. This provided 
a test of whether the RV candidates were in fact mostly 
face-on binaries (transiting planets prove that at least 
some RV candidate s are real planets — see, for example, 
Henry et al. (2000) — but this is not applic able to the 
vast majority ot systems). Han et al. (2001) concluded 
that most RV candidates did show orbital motions in 
the Hipparcos database, and w ere thus b i nary s tars, not 
planetary systems. However, Pourbaix (2001) showed 
the orbital analysis to be incorrect, and a proper statis- 
tical analysis reveals no credible detections; instead the 
results are consistent with randomly oriented orbits and 
most of the RV-detected objects being planetary in na- 
ture. 

A few RV-detected planetary systems have had their 
orbital geo metries constrained by Hubble Space Telescop e 
astrometry ( |Benedict et al |2002| [McArthur et al.|2010 ). 
Though impressive work, these do not represent discov- 
eries of new systems by astrometry, and the ratio of mea- 
surement precision to signal amplitude is low enough to 
make it unlikely the astrometry could have produced a 
detection by itself (or even have been made, given the 
time requirements of a blind search with Hubble) had 
the RV detection not already been present. |Pravdo et al.| 
( |2005[ ) successfully used astrometry to discover a brown 
dwarf companion to a M dwarf, a promising first step on 
the path to finding true planets. The brown dwarf was 



later confirme d by direct imaging dLloy d et al.||2006) 
Most recently, Pravdo fc Shaklanj (|2009[) claimed an as 



trometric detection ot a giant planet around a nearby M 
dwarf from the STEPS project, using standard CCD as- 
trometry from large aperture telescopes. However, this 
candidate was r apidly shown to be inconsistent with RV 
observations by Bean et al. (2010). 



It is thus with some trepidation that we announce the 
candidate substellar companions orbiting either the pri- 
mary or secondary stars in several binaries studied us- 
ing differential astrometry by PHASES — the Palomar 
High-precision Astrometric Search for Exoplanet Sys- 
tems. Given that other astrometrically "discovered" sub- 
stellar objects have not withstood the test of continued 
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observations, these may represent either the first such 
companions detected, or the latest in the tragic history 
of this challenging approach. 

Given the challenges of astrometry, why would it be 
considered a preferred way to detect planets in current 
and future searches? Astrometry has a number of ad- 
vantages over other techniques: 

1. Astrometry and RV provide information about the 
masses of companions to nearby stars. Other methods 
are insensitive to this fundamental property. 

• Because the reflex motion of the star is monitored, 
the mass of the companion is measured directly. For 
nearby systems that can be followed up by direct 
imaging, only RV and astrometry provide such in- 
formation. 

• The two-dimensional nature of the astrometric mea- 
surement provides unique mass estimates. In con- 
trast, RV detections only give the companion mass 
times the sine of the unknown inclination, M sini. 

2. Astrometry is effective in regimes where RV has re- 
duced precision: 

• Astrometry can operate over a wide range of stellar 
masses, rotational velocities, and spectral types, to 
better explore relationships between the properties 
of the host star and its planetary system. RV is 
most effective for slowly rotating, mid-to-late type 
stars. 

• Astrometric sensitivity increases with companion 
period, an opposite trend as RV. This is particularly 
important when identifying long period planets for 
direct imaging work, where the wider star-planet 
separation reduces technical challenges for imaging. 

• Astrometry is less se nsitive to surface vibra tions and 
star spots than RV ( |Makarov et al.||2009[ ). This is 
particularly important tor identifying the small (1 
/zas and 0.1 ms _1 ) signals of Earthlike planets in 
the habitable zones of nearby, Sunlike stars. This 
motivates future astrometric planet searches. 

• Astrometry is well suited to studying planets in bi- 
nary systems. This is the primary motivation for 
using the astrometric method for PHASES. RV can 
study planetary companions to a few binary sys- 
tems. For example, binaries with very large sky 
separations can be studied, which frequently implies 
large physical separations as well, and these evolve 
rather like single stars, revealing little new about 
planetary system formation and evolution. 

It is instead binaries with separations in the criti- 
cal ~ 10 — 50 AU range that can greatly contribute 
new information. This range is wide enough that 
planets can have stable orbits around either star if 
present, but close enough that the second star may 
influence formation of the planet in the first place. 
RV can study a few special cases of these binaries: 
those very close to the solar system ( e.g., 7 Cep; 



2010 1, and a few triple star systems, where a short- 



Campbell et al.|1988 Hatzes et al.|2003 ) so the com- 
ponents are spatially resolved, a tew high 



period stellar subsystem causes the spectral features 
to be split (e.g., the contro versial companion to 
HD 188753; jKonacki|[2005a[ ). In other cases, the 
stars are often both spatially unresolved and spec- 
trally blended, making precision velocities impossi- 
ble; even when the lines can be separated, precision 
RV on the double spectrum is challenging ( Konacki 



2005b). 



Thus, PHASES used astrometry to observe 51 binary 
systems with the goal of identifying new tertiary com- 
panions orbiting either of the bright stars in the system. 
Of these, 33 systems have more than 10 successful ob- 
servations, allowing a realistic chance for a companion 
search to be successful. Seven of those 33 systems are 
either triple or quadruple stars, five more are so distant 
that their physical separations fall well outside of the 50 
AU limit (two of these— HD 171779 and HD 221673— 
may have brown dwarf companions, and are presented 
in this paper), and six more have semimajor axes less 
than 10 AU, though these last two classes are useful 
to verify the astrometric technique. The remaining 15 
systems provide a sample from which the frequency of 
planets in closely separated binaries can be evaluated. 



Pfahl & Muterspaugh ( 2006 1 showed that stellar encoun- 
ters, even in star forming regions where the stellar den- 
sity is higher than typical space, are rare enough that 
only ~ 0.1% of closely separated binaries could pick up 
a planet that had not originally formed as a compan- 
ion in the binary itself but rather via an exchange or 
binary hardening event. Observed frequencies of plan- 
ets in close binaries that are higher than this value offer 
evidence of in situ formation. If giant planets do form 
in these binaries, it is likely the process must be rapid. 
Current core-accretion models predict slow formation, 
though the competing gravitational instability method 
shows promise at rapid formation. Thus, the frequency 
of planet formation in these close binaries evaluates the 
relative frequencies with which these (and other) modes 
of giant planet formation occur in nature. 

Unfortunately, the statistics of the number of binaries 
that have been observed by RV are difficult to evalu- 
ate. However, several planets have been found in close 
binaries by RV methods, certainly more than the 0.1% 
frequency predicted by non-in situ formation; see Table 
[l] The next challenge is to evaluate the planet frequency 
in a less biased manner. Though limited in size, the 
PHASES sample represents an attempt to contribute to 
this effort. 

This paper is the fifth in a series, analyzing the final re- 
sults of the PHASES project after its completion in late 
2008. The first paper describes the observing method, 
sources of measurement uncertainties, limits of observ- 
ing precisions, derives empirical scaling rules to account 
for noise sources beyond those predicted by the standard 
reduction algorithms, and presents the f ull catalog of 
astrometric m easurements from PHASES ( Muterspaugh 
et al.||2010d|). The second paper combines PHASES as- 



ponents are spatially resolved, a tew high contrast 
systems (such that the second star minimally im- 
pacts the spectrum; e.g., HD 126614; [Howard et al 



trometry with astrometric measurements made by other 
methods as well as RV observations (when available) to 
determine orbital solutions to the binaries' Keplerian mo- 
tions, determining physical properties such as component 
masses and system distance when possible ( |Muterspaugh| 
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Table 1 

Close Binaries with Substcllar Companions. 



System 



Object Type a o(AU) 



Mi/M 2 c R t (AU) d References 



HD 41004~ 
HD 12bbl4 

hp 188755*" 

HD iyb555~ 

III I 7 hi Ih I 

HD 221b'/3 



*Y Cephei 



18.5 



0.36 



GJ 8b 1 



1.59/0.34 



3.6 



1, 2 



rfD 410U4 



O.Y/1.0 



3, 4 ,b 



^0 



■F5?r 



-20 



U.V/U.4 



^5" 



0.4/0.7 



1.14b/0.324~ 



Y 



12.3 0.50 1.06/1.63" 



XT 



8, 9, 10 



1.3/U.b 



References. — (1) Campbell et al. 1988; (2) Hatzes et al. 2003; (3) Queloz et al. 2000; 
(4) Mugrauer & Neuhau ser 2005; (5) Lag range et al. 2006; (6) Zucker et al. 2004; (7) 
Howard et al. 2010; (8) |Konacki| j2005a^ (9) Eggenberger et al. 2007; (10) Mazeh et 
al. 2009; (11) Chauvin et al. 2006 



"p" indicates a giant planet companion and "bd" indicates the companion is a brown 



dwarf". 

t> WV.. 



When the eccentricity is unknown, the projected binary separation is used as an ap- 



rjraximataoj^except^ri^l^ 



and the binary itself has been resolved, leading to two possible solutions with a - 



10 



and 50 
T 



AU. 



cl The distance from the primary star at which a disk would be rapidly truncated by tides 
QPichardo et al.|20U5[>. 



; The companion star is a white dwart oi mass ~U.5AfQ. To estimate Ht at the time ot 
nrmaTToii7an*^ri^marroni^^ 



The com panion st ar itself is a binary w ith scmimajor axis .67 AU. Th is candidate is 
controversial due to minimal data in the discovery paper with sporadic observing cadence 
and a lack of evidence found by Eggenberger et al. (2007) and Mazeh et al. (2009). 



et al. [2010b I . The third paper presents limits on the ex- 
lstence of substellar tertiary companions, orbiting either 
the primary or secondary stars in those systems, that 
are found to be consistent with being simple binaries 
( Muterspaugh et al. 2010c). The fourth paper presents 
three- component orbital solutions to a known triple star 
system (63 Gem A = HD 58728) and a newly discovered 
triple sy stem (HR 2896 = HD 60318) |Muterspaugh et al. 
(|2010a|. Finally, the current paper presents candidate 



substellar companions to PHASES binaries as detected 
by astrometry. 

Astrometric measurements were made as part of the 
PHASES program at the Palom ar Testbed Interferom- 
eter (PTI " 
Palomar ]V 



Colavita et al. 19991, which was located on 
ountain near San Diego, California. It was 
developed by the Jet Propulsion Laboratory, California 
Institute of Technology for NASA, as a testbed for inter- 
ferometric techniques applicable to the Keck Interferom- 
eter and other missions such as the Space Interferometry 
Mission (SIM). It operated in the J (1.2/xm), H (1.6/zm), 
and K (2.2/im) bands, and combined starlight from two 
out of three available 40 cm apertures. The apertures 
formed a triangle with one 110 m and two 87 m baselines. 
PHASES observations began in 2002 continued through 
2008 November when PTI ceased routine operations. 

2. ALGORITHM FOR IDENTIFYING ASTROMETRIC 
COMPANIONS 

Blind searches were conducted to identify potential ter- 
tiary companions to_thePHA_SES binaries. An algorithm 
based on th at of Cumming et al.| ( |1999[ ) and Cumming 



et al. ( 2008 ) was modified tor use on astrometric data tor 



binary systems, as described in Paper III, and used to 
conduct blind searches for tertiary companions in these 
systems. 



The overall procedure is to create a periodogram of an 
F statistic comparing the goodness-of-fit x 2 between a 
single Keplerian model and that for a double Keplerian 
model for a number of possible orbital periods for the 
second orbit. The orbital periods selected were chosen 
to be more than Nyquist sampled, to ensure complete 
coverage, as P = 2fT/k where T is the span of PHASES 
observations, / = 3 is an oversampling factor, and k 
is a positive integer. Two searches were conducted for 
each binary: first, with the use of only the PHASES 
measurements, and second with both the PHASES and 
non-PHASES astrometry, to better constrain the wide 
binary motion during the search. In addition to the 
positive integer values of k, the period corresponding to 
k = 1/2 was evaluated to search for companions with 
orbits slightly longer than the PHASES span. 

The orbital period for which the F statistic peri- 
odogram has its maximum value is the most likely orbital 
period of a companion object. To ensure the peak is a 
real object rather than a statistical fluctuation, 1000 syn- 
thetic data sets with identical cadence and measurement 
uncertainties as the actual data were created and evalu- 
ated in the same manner. The fraction of these having a 
maximum F statistic larger than that of the actual data 
provided an estimate of the false alarm probability (FAP) 
that the signal is not caused by an actual companion. 

3. PHASES MEASUREMENTS 

PHASES differential astrometric measurements were 
obtained with the observing method and standard data 
analysis pipeline described in Paper I. The measurements 
themselves and associated measurement uncertainties are 
also tabulated in Paper I. The number of PHASES mea- 
surements available for each of the six systems being in- 
vestigated are listed in Table [2j 
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Table 2 

Number of PHASES and Non-PHASES Measurements and 
Unit Weight Uncertainties for Non-PHASES Measurements 



HD Number 


N P 


N P ,o 


N N p 


Nnp.o 


°>,° 




13872 


89 





103 


14 


0.013 


2.51 


171779 


54 





128 


12 


0.020 


2.79 


176051 


65 


1 


327 


12 


0.140 


4.48 


196524 


72 


1 


598 


48 


0.046 


2.78 


202444 


39 





286 


13 


0.123 


4.37 


221673 


98 


1 


333 


21 


0.056 


2.12 



Note. — The numbers of PHASES and non-PHASES 
astrometric measurements used for orbit fitting with each of 
the binaries being studied are presented in Columns 2 and 4 
respectively, along with the additional numbers of measure- 
ments rejected as outliers in Columns 3 and 5. Columns 6 
and 7 list the Itj measurement uncertainties for unit weight 
measurements from non-PHASES observations determined 
by iterating Keplerian fits to the measurements with removal 
of 3<r or greater outliers in either dimension. Columns 6 and 
7 are in units of arcseconds and degrees, respectively. 



4. NON-PHASES ASTROMETRY 

Measurements of binaries observed by PHASES made 
by previous astrometric techniques and ca taloged in the 
Washington Double Star Catalog (WDS Mason et al. 
2001 2010[) were assigned weights accor ding to the tor 



mula described by Hartkopf et al 
refined planet searches in which t 



pOOlj ). These allowed 
le binary orbit itself is 



better constrained by the longer duration, though lower 
precision, astrometric measurements. Including these 
measurements lifts some degeneracies in the double-orbit 
modeling. 

Unit weight uncertainties in separation and position 
angle were evaluated by the following iterative proce- 
dure. First guess values for the unit uncertainties of 24 
mas in separation and 1.8° in position angle were as- 
signed to the measurements of a given binary; these val- 
ues corresponded to previous experience u sing this pro- 
cedure on \x Ori ( Muterspaugh et aL]|2008 ). Second, the 
measurements were fit to a Keplerian model and the or- 
bital parameters were optimized to minimize the fit \ 2 - 
This intrinsically assumes the non-PHASES astrometric 
measurements are insensitive to the tertiary companions 
being sought, an assumption that will be justified given 
the small sizes of the perturbations detected. Third, the 
weighted scatter of the residuals in separation and posi- 
tion angle were evaluated, and the guessed unit uncer- 
tainties updated to make the rms scatter in each equal 
to unity. Fourth, the second and third steps were it- 
erated two more times, at which point the values con- 
verged. Fifth, the final unit uncertainties were multiplied 
by the square root of the reduced \ 2 (\/xf) 01 the fit, 
and refit one more time with these slightly larger weights. 
Sixth, if no residuals deviated by more than 3<r, the pro- 
cess ended, otherwise, the single measurement with the 
largest separation or position angle residual (weighted by 
its uncertainty) was flagged as an outlier, and removed 
from future fits. Seventh, the process was repeated at 
the first step. The resulting weights are listed in Tablep] 
and the measurements themselves are listed in Table |3] 



5. SUBSTELLAR COMPANIONS WITH HIGH DEGREES OF 
CONFIDENCE 

5.1. HD 176051 

HD 176051 (HR 7162, HIP 93017, WDS 18570+3254, 
and, though rarely used, the proper name of Inrakluk has 
been proposed) is an intriguing PHASES binary because 
its components are relatively low mass (1.07 and 0.71 
M©), and the system is relatively nearby ( 14.99 ± 0.13 
pc) a s determined by Hipparcos observations ( Soderhjelm 
1999| hereafter S99). Both of these qualities indicate as- 
trometric perturbations by tertiary c ompanions will have 
relativ ely large signals. The model of |Holman fc Wiegert 
(1999) (hereafter HW99) for determining which plane 



tary orbits in binary systems are stable long-term pre- 
dicts companions with periods as long as ~ 3000 days 
can have stable orbits. 

The initial PHASES-only search for companions found 
a most significant peak of z = 15.0 at a period of 581 days 
with FAP 0.0%. This low value inspired a revised search 
including both the PHASES and lower-precision non- 
PHASES astrometric observations to investigate whether 
the detection continued to be valid. The revised search 
finds the most significant peak of z = 63.1 at a period 
of 1004 days with FAP 0.0%. These two distinct orbital 
periods are probably aliases of each other and both are 
present in both periodograms; see Figure [T] 

The two peaks may indicate aliasing, orbital eccen- 
tricity, or confusion with the wide binary orbit. Both 
companion orbital periods were further explored using a 
double Keplerian model optimizing all orbital elements, 
including the companion orbital period and allowing for 
non-circular companion orbits. Both the fit % 2 and vi- 
sual inspection of the orbital solution confirmed that the 
longer period solution is more likely to be correct and 
the other is a harmonic. Additionally, smaller peaks cor- 
responding to 276, 225, and 7 days were also explored, 
but did not produce convincing solutions at all. 

The eccentricity of the subsystem orbit is not con- 
strained by the astrometry measurements. This is likely 
due to the PHASES measurements being high precision 
only in one dimension, making it difficult to use Kepler's 
second law to lift ambiguity between inclined circular or- 
bits and face-on eccentric ones. The eccentricity is fixed 
at zero in the present analysis. 

The best fit Keplerian stellar binary+circular subsys- 
tem orbital solution is presented in Table [4] and Figure [2j 
The substellar object is a planet 1.5 ±0.3 times the mass 
of Jupiter, assuming a distance of 15 pc and a stellar 
mass of 0.71 M , both based on the Hipparcos analysis 
by S99. If the planet is instead around the more massive 
star, the planet's mass would be twice as large. 

Interestingly, the binary and planetary orbits may be 
nearly coplanar. The mutual inclination $ of two orbits 
is given by 

cos 4> = cos i\ cos ?2 + sin i\ sin ii cos (fii — O2) (1) 

where i\ and 12 are the orbital inclinations and fii and 
£^2 are the longitudes of the ascending nodes. When RV 
measurements are not available, there exists ambiguity in 
which node is ascending, and two different values of the 
mutual inclination are possible (corresponding to Q\ — ^2 
varying by 180°). In this case, the two possibilities are 
$ = 18 ± 17 degrees or $ = 126.4 ± 6.6 degrees. 
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Table 3 

Non-PHASES Astrometric Measurements 



HD Number 


Date 


P 





Op 




Weight 


Outlier 




(year) 


(arcsec) 


(deg) 


(arcsec) 


(deg) 






13872 


1965.9100 


0.250 


49.70 


0.015 


2.81 


0.8 


1 


13872 


1966.1100 


0.230 


50.70 


0.029 


5.61 


0.2 





13872 


1966.7200 


0.220 


42.80 


0.012 


2.29 


1.2 





13872 


1967.0699 


0.190 


39.50 


0.016 


3.00 


0.7 





221673 


2006.7170 


0.560 


98.20 


0.102 


3.87 


0.3 





221673 


2006.9750 


0.554 


97.70 


0.030 


1.13 


3.5 





221673 


2007.2800 


0.580 


99.20 


0.040 


1.50 


2.0 





221673 


2008.8850 


0.540 


100.45 


0.018 


0.69 


9.5 






Note. — Non-PHASES astrometric measurements from the WDS Catalog are 
listed with la measurements uncertainties, and weights. Column 1 is the HD catalog 
number of the target star, Column 2 is the decimal year of the observation, Columns 
3 and 4 are the separation in arcseconds and position angle in degrees, respectively, 
Columns 5 and 6 are the 1-a uncertainties in the measured quantities from Columns 
3 and 4, Column 7 is the weight assigned to the measurement, and Column 8 is 1 if 
the measurement is a >3<r outlier and omitted from the fit, otherwise. (This table 
is available in its entirety in machine-readable and Virtual Observatory (VO) forms 
in the online journal. A portion is shown here for guidance regarding its form and 
content.) 
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Figure 1. Periodograms of the F statistic comparing models with and without tertiary companions to HD 176051 (HR 7162) in astrometric- 
only models. The left figure is for analysis only using the PHASES data, while the right is for combined analysis of PHASES and non- 
PHASES astrometric measurements. For HD 176051 the 1% FAP is at z = 8.37 for the PHASES-only analysis, and z = 11.33 for the 
combined analysis, as indicated by horizontal lines. 

the late summer/fall months of best weather at Palo- 
mar, and served as one of the easiest and most reliable 
PHASES targets to observe. As a result, 98 PHASES 
measurements were successfully taken of HD 221673. De- 
spite there being large amounts of data available, sin- 
gle Keplerian orbit fitting was frustrated from the early 
beginnings — the data show much more scatter than pre- 
dicted by the measurement uncertainties. 

The initial PHASES-only search for companions found 
a most significant peak of z = 32.9 at a period of 
1276 days (k = 9) with FAP 0.0%. This low value 
inspired a revised search including both the PHASES 
and lower-precision non-PHASES astrometric observa- 
tions to investigate whether the detection continued to 
be valid. The revised search finds a most significant peak 
of z — 95.9 at a period of 1435 days (k — 8, within one 
sample of the peak value for the PHASES-only search) 



5.2. HD 221673 

HD 221673 (72 Peg, HR 8943, HIP 1 16310 
23340+3120) is a pair of mid K giants. |Baize| ( [1962 1 
flagged it as possibly containing a variable star based on 
the scatter in the differential magnitude measurements 
by various observers. However, Hipparcos photometry 
shows a scatter of only 6 mmag. T he revised Hipparco s 



based parallax is 5.94 ± 0.45 mas ( |van Leeuwen||2007[ ). 
This parallax and the best fit single Keplerian model 
predict an average stellar mass of 2 Mq . Its long orbital 
period (~ 800 years) implies a large range of orbits in 
which companions can be stable — up to 47 years accord- 
ing to the criteria of HW99. Thus, binary dynamics are 
expected to have a smaller impact on planet formation 
in this system than others. 

HD 221673 is extremely bright at infrared wavelengths 
(K = 1.76), is observable for long stretches during 



G 



Table 4 

Orbit Model for HD 176051 



Table 5 

Orbit Model for HD 221673 



Parameter 


Value 


Uncertainty 


Parameter 


Value 


Uncertainty 


Pa d (Hpjvs) 

1 A— a v uo, J a / 


22430 


15 


P a d ( davsl 

J\ — B V J / 


179811 


27745 




41384 


23 


Ta h ('MHJD') 


16818 


3658 




0.2667 


0.0022 


B 


0.322 


0.047 


a,A-B (arcsec) 


1.2756 


0.0023 


cia-b (arcscconds) 
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41 
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4.0 


^Ba-Bb (deg) 





(Fixed) 


^Ba-Bb (deg) 





(Fixed) 
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11 
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Note. — Best fit orbital elements in the 
Campbell basis for HD 176051, with la uncer- 
tainties. 
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Figure 2. Time series of PHASES observations of HD 176051 
(HR 7162) for the 1016 day subsystem orbital solution, measured 
along an axis at angle 167° from increasing differential right ascen- 
sion through increasing differential declination; it was along this 
axis that the PHASES measurements are typically most sensitive. 
For clarity, only measurements with uncertainties along this axis 
of 100 /tas along this axis are shown. 

with FAP 0.0%. The periodograms are plotted in Figure 

El 

A few single-compo nent RV measurements of 72 Peg 
have been publish ed by Tokovinin & Smekhov ( 2002[ ) and 
Abt et al. (19801. However, the relatively small number 
and short time coverage of each RV data set reduces any 
impact they have on constraining either the binary or- 
bit or confirming the existence of additional components 
(especially low mass, long period companions). Neither 
set shows variation in the velocities and as a result is not 
used in the orbital analysis. 

A refined fit to the astrometric measurements was 
made for each of the three most significant peaks in the 
periodogram — the companion orbital period was seeded 
with values of 1435, 478, and 205 days. Both circular and 
full Keplerian models were attempted as the companion 
orbit for each of the three periods being explored. The 



Note. — Best fit orbital elements in the 
Campbell basis for HD 221673, with lcr un- 
certainties. 

longest period corresponded to the best fit of the three, 
though the companion's eccentricity was not constrained. 
The best fit circular model converged with an orbital pe- 
riod of 1539 days and \ 2 = 1549.2, with 850 degrees of 
freedom, and is presented in Tableland Figure [3j While 
X 2 is larger than the number of degrees of freedom, it is 
significantly improved compared to the model without a 
tertiary companion, for which \ 2 = 2669.9 with 855 de- 
grees of freedom. The 1539 day companion is a brown 
dwarf with 35 ± 4 times the mass of Jupiter. 

However, the remaining scatter and presence of other 
peaks in the periodogram (especially that at 478 days) 
leads one to question whether adding yet another Ke- 
plerian representing a fourth component to the system 
would yet further improve the fit. A 3-Keplerian fit 
was seeded with the best parameters from the three- 
component, 2-Keplerian model, as well as the best or- 
bit for a 478 day companion as determined by the ini- 
tial search. First, both subsystem orbits were assumed 
to be circular. This led to an improvement in the fit 
from x 2 = 1549.2 with 850 de grees of freedom for the 
2-Keplerian model to \ 2 = 1422.9 with 845 degrees of 
freedom for that with three Keplcrians summed by su- 
perposition. This is only a modest improvement, and 
at this time no detection is claimed for a second unseen 
object. Alternatively, the remaining scatter could be re- 
lated to the 6 mmag photometric variability. 

6. SUBSTELLAR COMPANIONS WITH REDUCED LEVELS 
OF CONFIDENCE 

6.1. HD 13872 

HD 13872 (21 Ari, HR 657, HIP 10535, WDS 
02157+2503) is a bright star with mid-F dwarf spectrum. 
In |1967[ it was realized to be a visual binary system with 
separation les s than an arcsecond and roughly equal lu- 
minosities by Couteau Since its first orbit determina- 
tion, there have been questions as to whether its spec- 
tral ty pe fit the total system m ass as measured by the 
orbit. |Couteau fc Morel (1982) proposed that an ad- 
ditional component must exist in the system to explain 
the overly large total mass. However, their estimate for 
the total mass was in err or due to too sma ll a value of 
the parallax (of 15 mas). Tokovinin (1987) pointed out 
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Figure 3. Pcriodograms of the F statistic comparing models with and without tertiary companions to HD 221673 (72 Peg) in astrometric- 
only models. The left figure is for analysis only using the PHASES data, while the right is for combined analysis of PHASES and non- 
PHASES astrometric measurements. For HD 221673 the 1% FAP is at z = 7.56 for the PHASES-only analysis, and z = 13.96 for the 
combined analysis, as indicated by horizontal lines. 
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plotted in Figure [5] 
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Figure 4. Time series of PHASES observations of HD 221673 
(72 Peg) for the 1539 day subsystem orbital solution, measured 
along an axis at angle 155° from increasing differential right ascen- 
sion through increasing differential declination; it was along this 
axis that the PHASES measurements were typically most sensitive 
for this system. For clarity, only measurements with uncertainties 
along this axis of 100 fi&s or less are shown. 

that a parallax of 20.1 mas would give a normal sum of 
masses, a p arallax later confirmed by Hipparcos ( |Perry- 
man et al. 19971. However, this is only the case if the 
star's mass ratio is near unity, which had not been deter- 
mined previous to the current investigation, in which a 
mass ratio Mb /Ma — 1.027 ± 0.032 has been measured. 

Fifty-one RV measurements of HD 13872 were made 
with Tennessee State Unive rsity's 2 m Automated Spec- 
troscopic Telescope (AST; |Eaton k. Williamson||2007[ ) 
and echelle spectrograph to obtain the stellar mass ratio 
and better constrain the binary orbit. These measure- 
ments are listed in Table [6l The standard data reduction 
pipeline for determining binary star velocities from AST 
data was used, as described in Paper II. The RV orbit is 



Figure 5. The RV orbit of the 21 Ari binary with measurements 
from TSU's AST. These measurements enable the binary mass ra- 
tio and component masses to be measured, but lack the precision 
necessary to confirm the presence of a planetary companion. The 
velocities of both components are shown in the top graph, where 
the individual components A and B are labeled. The middle graph 
shows the measurement residuals for component A, and the bottom 
graph those for component B. 

The initial PHASES-only search for companions found 
a most significant peak of z = 6.73 at a period of ~ 770 
days with FAP 3.0%. This relatively low value inspired a 
revised search including both the PHASES and lower pre- 
cision non-PHASES astrometric observations, to investi- 
gate whether the detection continued to be valid. The 
revised search finds a most significant peak of z = 7.49 at 
a period of - 1284 days with FAP 0.3%. The non-zero 
FAP, especially when only the PHASES measurements 
are considered, prevents identification of this as being 
counted among the strongest of candidates, but is an in- 
triguing possibility, since it would correspond to a giant 
planet. The periodograms are plotted in Figure [6j 

Orbit fitting was refined by allowing the companion 
orbital period to be optimized, non-circular orbits to be 
considered, and by adding to the data set the 51 two- 
component RV measurements from the TSU AST spectra 
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Figure 6. Periodograms of the F statistic comparing models with and without tertiary companions to HD 13872 (21 Ari) in astrometric- 
only models. The left figure is for analysis of only the PHASES data, while the right is for combined analysis of PHASES and non-PHASES 
astrometric measurements. For HD 13872 the 1% FAP is at z = 7.33 for the PHASES-only analysis, and z = 6.27 for the combined analysis, 
as indicated by horizontal lines. 

that span 358 days, enabling a full three-dimensional or- 
bit to be evaluated, including the distance to the system 
and component masses. There is a small, but insignifi- 
cant, fit improvement if the companion is assumed to be 
around star A instead of star B. The eccentricity of the 
subsystem orbit is not well constrained due to the small 
signal size and the majority of the highest precision as- 
trometry being along only one axis on the sky. Thus, the 
eccentricity is fixed at zero in the present analysis. The 
2-Keplerian orbit model is presented in Table [7] If real, 
the perturbation corresponds to a giant planet of mass 
1.40 ± 0.36 Mj. The reflex motion of star A due to the 
presence of the companion is plotted in Figure [7] with the 
A-B binary orbit removed. 

The model of HW99 predicts only planets with orbital 
periods less than 210 days will be stable. This is far 
smaller than the value of 925 ± 90 days that best fits the 
combined measurements. While the model from HW99 
is broad in scope for generalized orbits, an analysis spe- 
cific to the configuration of HD 13872 would be beneficial 
to explore whether there are additional islands of stabil- 
ity for companion orbital periods. However, the mutual 
inclination of the binary and planet orbits is not con- 
strained well, limiting the utility of system-specific sta- 
bility analysis. Clearly this cannot be considered a high 
confidence detection at this time. 

Because the data are not of high enough quality to 
constrain the companion eccentricity, the FAP is greater 
than 0.1%, and the planet's orbital period may be unsta- 
ble according to the criteria of HW99, the reality of this 
planet is highly uncertain. HD 13872 will be an interest- 
ing object for continued study, but in this case it might 
not be surprising if future observations do not confirm 
the presence of a giant planet in the system. 

6.2. HD 202444 

HD 202444 (r Cyg, 65 Cyg, HR 8130, HIP 104887, 
WDS 21148+3803) is classified as an early F subgiant 
with a G dwarf companion. Various reports have sug- 



gested that the primary is a S Scuti or 7 Doradus vari- 
able, though these have not been confirmed. 

The search for planetary companions to HD 202444 
is more complicated than for other stars presented in 
this paper. It appears that a companion object may 
exist with an orbital period comparable to the span of 
PHASES observations (only two observations were taken 
outside of the 1155 day span from MJD 53234-54389 
when most observations of reasonable cadence were made 
whereas the companion orbital period is over 800 days). 
Because the companion search software reoptimizes both 
the wide binary orbit model and the perturbing model 
every time a fit is made, the signal could be absorbed into 
that of the wider binary when only the shorter timespan 
PHASES data were analyzed. Thus, no compelling evi- 
dence for a companion was present when only PHASES 
measurements were analyzed — the initial PHASES-only 
search for companions found a most significant peak of 
z = 5.93 at a period of 25.5 days with FAP 19.1%. Analy- 
sis of the combined PHASES and non-PHASES data sets 
showed a larger value of x 2 than one would have antici- 
pated based on fits to the individual data sets, prompt- 
ing a second search for tertiary companions, this time 
using all the astrometric measurements. When the non- 
PHASES measurements were added to the analysis, the 
extended coverage of the binary orbit prevented much 
of the ability to adjust the binary orbit to include the 
perturbations caused by possible companions. The re- 
vised search finds a very significant peak of z = 51.85 at 
a period of 826 days with FAP 0.0%. The orbit stabil- 
ity criteria of HW99 predict that companions with pe- 
riods shorter than 2200 days are stable in HD 202444, 
which includes all candidate periods identified in the pe- 
riodograms. The periodograms are plotted in Figure [51 

The orbit fitting was refined by seeding a full Keple- 
rian fit with the best orbital parameters corresponding 
to the three largest peaks in the full data periodogram at 
periods 826, 534, and 252 days. The best fit occurred for 
the longest of these periods, for which the eccentricity of 
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Figure 7. Motion of the center of light of the 925 day subsystem in HD 13872 (21 Ari) measured along an axis at angle 156°, measured 
from increasing differential right ascension through increasing differential declination; it was along this axis that the PHASES measurements 
were typically most sensitive for this binary. For clarity, only measurements with uncertainties along this axis of 100 /xas or less are shown. 
Left: the center-of-light motion as a function of Modified Julian Date. Right: the center-of-light motion phase-wrapped about the orbital 
period and plotted covering two cycles. 
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Figure 8. Periodograms of the F statistic comparing models with and without tertiary companions to HD 202444 (r Ceti) in astrometric- 
only models. The left figure is for analysis only using the PHASES data, while the right is for combined analysis of PHASES and 
non-PHASES astrometric measurements. For HD 202444 the 1% FAP is at z = 8.92 for the PHASES-only analysis, and z = 10.12 for the 
combined analysis, as indicated by horizontal lines. 



the subsystem orbit was constrained to be 0.43 ± 0.17, 
so the full Keplerian model is accepted. The final fit has 
a subsystem period of 810 days, and the fit x 2 is 744.7 
with 636 degrees of freedom, compared to 994.7 and 643, 
respectively, for the single Keplerian fit. The best-fit 
two-Keplerian model parameters are presented in Table 
1 _ 

Since the companion is not detected when only the 
PHASES measurements or non-PHASES measurements 
are analyzed individually, there are reasons to doubt 
the authenticity of this proposed companion. Evidence 
for the companion only appears when the high-precision 
measurements are coupled with the measurements span- 
ning more time. The reflex motion of the star in the 



subsystem is plotted in Figure [9j in which four measure- 
ments with uncertainties larger than 100 fias projected 
onto the selected axis are not shown; most of these points 
are consistent with this fit, though one suppressed mea- 
surement at MJD 53285 is a 5.4er outlier. This mea- 
surement was taken with PTFs less reliable (and infre- 
quently used) south-west baseline, rather than the stan- 
dard north-south baseline. If the companion is a real 
object, it is at the border between the realm of brown 
dwarfs and giant planets, with a mass of 12.3 ± 2.3 Mj, 
for which a stellar mass of 1.36 Af Q and a distance of 
20.37 ± 0.25 pc have been assumed based on the work 
of S99. In case the candidate is not real, the single Ke- 
plerian binary-only orbit model was presented in Paper 
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Table 6 

AST Velocities of 21 Ari 



Table 7 

Orbit Model for HD 13872 



Day 


RV A 


&RV A 


RV B 


a RV,B 


(HMJD) 


(kms -1 ) 


(kms Ll ) 


(kms^ 1 ) 


(kins" 1 ) 


54884.158 


-56.07 


0.31 


-37.13 





50 


54887.098 


-55.91 


0.31 


-37.31 





50 


54888.089 


-56.39 


0.31 


-37.68 





50 


54890.089 


-56.22 


0.31 


-37.55 





50 


54891.094 


-56.63 


0.31 


-37.83 





50 


54892.111 


-56.05 


0.31 


-37.57 





50 


54893.096 


-56.27 


0.31 


-37.96 





50 


54898.123 


-56.45 


0.31 


-38.16 





50 


54902.111 


-55.84 


0.31 


-37.67 





50 


54903.111 


-56.02 


0.31 


-38.15 





50 


54904.125 


-56.40 


0.31 


-38.06 





50 


54906.110 


-56.26 


0.31 


-38.58 





50 


54908.110 


-56.11 


0.31 


-37.83 





50 


54909.117 


-56.48 


0.31 


-38.54 





50 


54910.110 


-56.00 


0.31 


-37.81 





50 


54982.467 


-54.77 


0.31 


-39.88 





50 


54995.464 


-54.21 


0.31 


-40.33 





50 


55008.401 


-53.80 


0.31 


-40.84 





50 


55021.463 


-53.08 


0.31 


-40.59 





50 


55032.425 


-52.77 


0.31 


-41.16 





50 


55045.389 


-52.61 


0.31 


-41.96 





50 


55052.287 


-51.60 


0.31 


-41.63 





50 


55061.360 


-51.35 


0.31 


-42.48 





50 


55066.468 


-51.19 


0.31 


-42.59 





50 


55080.285 


-50.89 


0.31 


-42.69 





50 


55083.447 


-50.71 


0.31 


-42.61 





50 


55092.313 


-49.87 


0.31 


-43.20 





50 


55093.263 


-50.51 


0.31 


-42.75 





50 


55094.434 


-50.26 


0.31 


-43.87 





50 


55096.357 


-49.96 


0.31 


-43.80 





50 


55099.413 


-49.60 


0.31 


-44.33 





50 


55104.390 


-49.10 


0.31 


-44.73 





50 


55105.390 


-49.06 


0.31 


-44.09 





50 


55106.390 


-48.42 


0.31 


-45.48 





50 


55113.390 


-48.31 


0.31 


-45.75 





50 


55118.365 


-48.90 


0.31 


-45.09 





50 


55119.365 


-48.50 


0.31 


-45.18 





50 


55120.365 


-48.25 


0.31 


-45.56 





50 


55124.340 


-48.17 


0.31 


-45.34 





50 


55126.340 


-48.46 


0.31 


-45.73 





50 


55137.315 


-47.84 


0.31 


-46.13 





50 


55139.440 


-47.00 


0.31 


-46.51 





50 


55141.315 


-47.21 


0.31 


-46.91 





50 


55145.290 


-47.08 


0.31 


-46.75 





50 


55153.258 


-47.00 


0.31 


-47.01 





50 


55158.265 


-46.91 


0.31 


-46.63 





50 


55161.240 


-46.84 


0.31 


-46.38 





50 


55241.131 


-42.89 


0.31 


-51.20 





50 


55241.155 


-42.63 


0.31 


-52.18 





50 


55242.131 


-42.60 


0.31 


-51.61 





50 


55242.140 


-42.73 


0.31 


-52.00 





50 


Note. - 


Two-component RV measurements of 21 Ari 



Parameter 



Value Uncertainty 



(HD 13872) from the AST. 

II. Continued high-precision observations of the system 
spanning at least five years would help clarify the situa- 
tion. 

7. SUBSTELLAR COMPANION CANDIDATES WITH 
AMBIGUOUS OR UNCERTAIN ORBITAL 
CHARACTERISTICS 

7.1. HD 171779 

HD 171779 (HR 6983, HIP 91013, WDS 18339+5221) 
was the first binary for which PHASES observati ons were 
published as a dem onstration of the technique (Lane & 
Muterspaugh 2004). In total, 54 differential astrometric 
measurements were made at PTI of this pair of K giant 
stars. 

The initial PHASES-only search for companions found 



Pa-b (days) 


8622.7 


4.4 


T A _fl (MHJD) 


46497.2 


4.3 


e A-B 


0.68119 


0.00096 


M Aa+Ab (M ) 


1.338 


0.032 


M B (Mq) 


1.374 


0.027 


*A-b (deg) 


104.437 


0.025 


uiA-B (deg) 


263.927 


0.031 


^A-B (deg) 


55.823 


0.032 


P Aa -Ab (days) 


925 


90 


T Aa _ Ab (MHJD) 


54092 


62 


6Aa-Ab 





(Fixed) 


M Ab /M Aa 


0.00100 


0.00023 


L A b/L Aa 





(Fixed) 


*Aa-Ab (deg) 


71 


45 


^Aa-Ab (deg) 





(Fixed) 


f^Aa-Ab.i (deg) 


211 


55 


Vb AST (km s _1 ) 


-46.892 


0.053 


d (pc) 


48.90 


0.33 


X 2 and dof 


406.1 


471 



Note. — Best-fit orbital elements in the 
Campbell basis for HD 13872, with la un- 
certainties. 



Table 8 

Orbit model for HD 202444 



Parameter 


Value 


Uncertainty 


Pa-b (days) 


18125.4 


7.7 


T A _ B (MHJD) 


47553 


17 


SA-B 


0.2392 


0.0012 


a A _g (arcseconds) 


0.9130 


0.0013 


i A -B (deg) 


134.44 


0.15 


ui A ~B (deg) 


298.77 


0.19 


(deg) 


339.75 


0.13 


-PBa-Bb (days) 


810 


18 


T Ba -Bb (MHJD) 


53139 


48 


EBa-Bb 


0.43 


0.17 


acOL (/ias) 


796 


119 


*Ba-Bb (deg) 


92.6 


1.9 


^Ba-Bb (deg) 


90 


19 


^Ba-Bb,l (deg) 


78.7 


2.5 


X 2 and dof 


744.7 


636 



Note. — Best-fit orbital elements in the 
Campbell basis for HD 202444, with la un- 
certainties. 

a most significant peak of z = 7.76 at a period of 
1663 days with FAP 1.4%. This low value inspired a 
revised search including both the PHASES and lower- 
precision non-PHASES astrometric observations to in- 
vestigate whether the detection continued to be valid. 
The revised search finds a most significant peak of z = 
18.85 at a period of 2328 days with FAP 0.0%. The pe- 
riodogra ms r esulting from these searches are presented 
in Figure |10| The two approaches find orbital periods 
within two samples of each other at k = 5 and k = 7, 
within the / = 3 oversampling factor of the search peri- 
ods. Thus, the same signal is detected with both ap- 
proaches. The stability criteria established by HW99 
predict orbital periods up to 5200 days or more would 
be stable in this system. The best-fit Keplerian+circular 
subsystem orbital solution is presented in Table [9] and 
Figure [TTj 

However, this is not the only significant peak in the 
periodogram. In an attempt to uncover which of these 
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Figure 9. Motion of the center of light of the 810 day subsystem in HD 202444 (r Cyg) measured along an axis at angle 171° degrees, 
measured from increasing differential right ascension through increasing differential declination; it was along this axis that the PHASES 
measurements were typically most sensitive. For clarity, only measurements with uncertainties along this axis of 100 /xas or less are shown. 
Left: the center-of-light motion as a function of Modified Julian Date. Right: the center-of-light motion phase-wrapped about the orbital 
period and plotted covering two cycles for continuity. 
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Figure 10. Pcriodograms of the F statistic comparing models with and without tertiary companions to HD 171779 (HR 6983) in 
astrometric-only models. The left figure is for analysis of only the PHASES data, while the right is for combined analysis of PHASES and 
non-PHASES astrometric measurements. For HD 171779 the 1% FAP is at z = 7.90 for the PHASES-only analysis, and z = 8.36 for the 
combined analysis, as indicated by horizontal lines. 



might be a true signal, further efforts were made to refine 
the orbital fits with the subsystem orbital period seeded 
near 2328 days (k = 5), 647 days (k = 18), 448 days 
(k = 26), 831 days (k = 14), and 233 days (k = 50). 
In all of these cases except for the 448 day selection, 
the eccentricity of the subsystem orbit could not be con- 
strained by the astrometric measurements; in the case 
of the 448 day period, the best-fit non-zero eccentricity 
was 0.91 ± 0.34, which is only poorly defined. Further- 
more, visual inspection of the orbital fit to all but the 
2328 day perturbation shows the others are almost cer- 
tainly a result of observin g ca dence rather than a true 
periodic signal — see Figure 12 The 2328 day signal (re- 



fined to 2324 days after full orbit fitting) is slightly longer 
than the 1940 day span of PHASES observations, and the 
majority of the measurements might be best represented 
as a linear trend, though a circular orbital solution is 
also possible. If this corresponds to a real companion, 
the substellar object is either a very massive planet or 
a brown dwarf roughly 10 times the mass of Jupiter, as- 
suming a distance of 196 pc (based on the revised Hip- 
parcos parallax) and a stellar mass of 1.4 Mq (derived 
from the binary orbit). 

The multiple peaks in the periodogram, potential for 
aliasing, inability to constrain the eccentricity of the 
companion, and best-fit orbit period being longer than 
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Table 9 

Orbit model for HD 171779 



Parameter 


Value 


Uncertainty 


Pa d 1 HrlVR 1 


75200 


2464 




21156 


295 


&A — B 


0.4161 


0.0083 


a,A-B (arcscconds) 


0.2524 


0.0072 


iA-B (deg) 


48.0 


1.4 


^A-B (deg) 


262.4 


2.9 


flA-B (deg) 


57.5 


1.3 


^Ba-Bb (days) 


2324 


250 


T Ba -Bb (MHJD) 


53375 


21 


eBa-Bb 





(Fixed) 


a-COL (/^as) 


160 


48 


*Ba-Bb (deg) 


66 


21 


w Ba -Bb (deg) 





(Fixed) 


^Ba-Bb,i (deg) 


157 


33 


X 2 and dof 


335.8 


352 



Note. — Best-fit orbital elements in the 
Campbell basis for HD 171779, with la un- 
certainties. 

the PHASES observation span prevent this detection 
from having high confidence. This system warrants 
further investigation over longer timespans to evaluate 
whether the detected perturbation is real and to con- 
strain the orbit. 




52925 53290 53655 54020 54365 54750 
t(MJD) 

Figure 11. Time series of PHASES observations of HD 171779 
(HR 6983) for the 2328 day subsystem orbital solution, measured 
along an axis at angle 175° from increasing differential right ascen- 
sion through increasing differential declination; it was along this 
axis that the PHASES measurements are typically most sensitive. 
For clarity, only measurements with uncertainties along this axis 
of 100 fi&s or less are shown. 



7.2. HD 196524 

HD 196524 (/? Del, 6 Del, HR 7882, HIP 101769, WDS 
20375+1436) is a pair of F5 subgiants. The system is the 
brightest star in its constellation, despite being given the 
Bayer designation /3. It and the fainter a Del were given 
proper names in the mid 1800 's by Niccolo Cacciatore 
when he compiled the Palermo Star Catalogue; a and /? 
Del have since been cataloged as Sualocin and Rotanev, 
respectively. These names are peculiar because they are 
the reverses of Nicolaus and Venator, the Latinized ver- 



sions of Cacciatore's own names 

The initial PHASES-only search tor companions found 
a most significant peak of z = 7.76 at a period of 6.81 
days with FAP 1.9%. This value is somewhat suspi- 
cious since PHASES observations were often scheduled 
the same nights each week. However, a nearly equal 
height peak of z = 7.53 occurs at period 422 days, and 
yet another with z = 7.46 at 203 days. The presence 
of three peaks at very different orbital periods, potential 
for aliasing confusion in the peaks, and low FAP values 
inspired a revised search, including both the PHASES 
and lower-precision non-PHASES astrometric observa- 
tions, to explore the sample of companion orbital peri- 
ods to investigate whether the detection continued to be 
valid. The revised search finds a most significant peak of 
z = 16.77 at a period of 439 days with FAP 0.0%, within 
one sampling of the secondary peak in the original anal- 
ysis. Furthermore, the peaks at ~ 200 and 6.81 days are 
still present, with values of z = 15.88 and z = 12.88, 



respectively. The periodograms are plotted in Figure 13 
Fits to double Keplerian models were seeded at the 
three potential companion orbital periods to explore 
which converged on the more satisfactory solution. The 
longest period model converged to a final orbital period 
of 435 days for a circular model with % 2 = 1551.4 and 
1328 degrees of freedom. The full Keplerian version was 
unable to constrain the eccentricity and the fit failed. 
Analysis of the middle period solution was able to con- 
strain modestly the eccentricity to e = 0.67 ±0.24 with a 
final x 2 = 1549.1 and 1326 degrees of freedom. Finally, 
the solution with a period near one week found a best-fit 
period of 6.8116 ± 0.0016 days for both circular and ec- 
centric (e = 0.27 ± 0.28) models, with X 2 = 1570.5 and 
X 2 = 1570.2, respectively. 

Though 96 spectra of this system have been obtained 
by TSU's AST, the spectral features of the two compo- 
nents were blended in all cases and could not be used 
for additional analysis. However, this does make it less 
likely that the 6.81 day signal is evidence of a real com- 
panion. Th e two longe r p erio d so lutions are presented 

The Hipparcos based 



in Table [10] and Figures [14] and 15 
parallax of 32.5 ± 0.7 mas and average component mass 
of 1.67 Mq from S99 can be used to convert the stellar 
reflex motion to companion mass. If real, the 435 day 
companion would be a giant planet of 9 ± 1.6 times the 
mass of Jupiter whereas the 202 day companion would 
be a brown dwarf of 15 ± 3.9 times the mass of Jupiter. 

At this point, it is not possible to distinguish whether 
either solution represents an actual companion, nor 
which model is preferred. Adding to the challenge of 
evaluating this system is the finding that for both so- 
lutions the perturbation orbit has an orientation on the 
sky that is more closely aligned with the major axis of 
the typical PHASES uncertainty ellipse. Where the sig- 
nal is strongest, the astrometric precision is the worst. 
HD 196524 is a system that would benefit from contin- 
ued observation by future astrometric efforts, especially 
those capable of truly two-dimensional measurements or 
that are sensitive to the perpendicular axis compared to 
PHASES. 



8. CONTINUED STUDIES 

Candidate substellar objects discovered by PHASES 
astrometry include: 
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Figure 12. Time series of PHASES observations of HD 171779 (HR 6983) for two candidate orbital periods for companions to HD 171779, 
showing these solutions are likely due to cadence rather than an actual object. The best-fit model for the wide A-B binary motion based on 
a simultaneous fit with the perturbing model has been removed in each case, leaving only the motion of the center of light of the subsystem. 
Measurements are plotted along an axis at angle 175°, measured from increasing differential right ascension through increasing differential 
declination; it was along this axis that the PHASES measurements were most sensitive for this binary. For clarity, only measurements 
with uncertainties along this axis of 100 ^tas or less are shown. Left: the candidate orbital period is 628 days. Right: the candidate orbital 
period is 451 days. 
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Figure 13. Pcriodograms of the F statistic comparing models with and without tertiary companions to HD 196524 (Rotancv//3 Del) in 
astrometric-only models. The left figure is for analysis of only the PHASES data, while the right is for combined analysis of PHASES and 
non-PHASES astrometric measurements. For HD 196524 the 1% FAP is at z = 8.37 for the PHASES-only analysis, and z = 8.76 for the 
combined analysis, as indicated by horizontal lines. 



a planet slightly more massive than Jupiter around one 
of the stars in HD 176051, 

one or more brown dwarfs around HD 221673, 

a possible Jovian planet orbiting one of the stars in the 
HD 13872 system, though this has low confidence given 
the non-zero FAP of the signal and the prediction that 
the orbit may not be stable over long periods of time, 

a possible very massive planet or low mass brown dwarf 
orbiting one of the stars in the HD 202444 binary, 
though this detection has reduced confidence because 



detection is not possible from PHASES measurements 
alone, but only reveals itself when lower precision as- 
trometry covering longer time periods aid in constrain- 
ing the binary orbit, 

a possible very massive planet or brown dwarf compan- 
ion in the HD 171779 system, though at present it is 
impossible to distinguish which of several possible or- 
bital periods are correct, and 

a possible massive planet or low-mass brown dwarf 
companion orbiting one of the stars in the HD 196524 



14 




Orbital Phase 



1 

o.e 

0.6 
0.4 

¥ - 2 

£ 
O-0.2 
-0.4 
-0.6 
-0.6 
-1 



i 


1 


1 


■ 












. J 




i 




K 


. r '- - - 










J? 

























0.5 



1 

Orbital Phase 



1.5 



Figure 14. Phase-wrapped PHASES observations of HD 196524 (Rotanev//3 Del) for the longest candidate orbital period for a tertiary 
companion. The best-fit model for the wide A-B binary motion based on a simultaneous fit with the perturbing model has been removed, 
leaving only the motion of the center of light of the subsystem. Left: motion along the right ascension axis; for clarity, only measurements 
with uncertainties projected on the right ascension axis 200 fias or smaller are shown. Right: motion along the declination axis; for clarity, 
only measurements with uncertainties projected on the declination axis 100 /jas or smaller are shown. 
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Figure 15. Phase-wrapped PHASES observations of HD 196524 (Rotanev//3 Del) for the 202 day candidate orbital period for a tertiary 
companion. The best-fit model for the wide A-B binary motion based on a simultaneous fit with the perturbing model has been removed, 
leaving only the motion of the center of light of the subsystem. Left: motion along the right ascension axis; for clarity, only measurements 
with uncertainties projected on the right ascension axis 200 ^tas or smaller are shown. Right: motion along the declination axis; for clarity, 
only measurements with uncertainties projected on the declination axis 100 /jas or smaller are shown. 



system, though which of at least two possible orbital 
periods are correct cannot be determined at this time. 

Of the 15 binary PHASES targets observed 10 or more 
times having semimajor axis in the 10-50 AU range, Pa- 
per III demonstrates the present data set cannot rule out 
planetary mass companions in any stable orbit for 4. Of 
the remaining 11 systems, there is strong evidence for a 
Jovian planet companion to HD 176051, while HD 13872 
may also host a planet, though this detection is with 
lower confidence. Furthermore, the remaining nine sys- 
tems in which no companions were detected but for which 
the constraints included some planetary mass compan- 
ions (HD numbers 5286, 76943, 81858, 114378, 137107, 



140436, 202444, 207652, and 214850) have a large range 
of unexplored orbital periods for which giant planetary 
companions cannot be ruled out. This implies that ei- 
ther the PHASES program was incredibly lucky, or giant 
planets are fairly common in close binary systems. The 
growing number of such systems being detected that are 
listed in Table [l] suggests that the latter explanation is 
more likely. 

|Thebault et al.| ( |2004 ) examined the formation of 7 
Cephei's gas g iant planet in the core-accretion scenario 



( Mizuno[ 
tions of t 



1980), subject to the gravitational perturba- 
le binary companion on a moderately eccentric 



(e = 0.36) orbit. Assuming a massive gaseous disk, they 
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Table 10 

Orbit Models for HD 196524 



Parameter 


Value 


Uncertainty 


Value 


Uncertainty 


r A — B V uc v»j 


9745.6 


1.4 


9745.8 


1.4 




37960 

•-J 1 <J\J \J ■ \J 


4.1 


37961 5 


4.0 


^A — B 


0.35632 


0.00070 


0.35595 


0.00069 


cia-b (arcscconds) 


0.43676 


0.00016 


0.43701 


0.00016 


iA-B (deg) 


61.289 


0.035 


61.323 


0.030 


WA-s (deg) 


168.81 


0.14 


168.86 


0.13 


&A-B (deg) 


357.206 


0.033 


357.179 


0.029 


-PBa-Bb (days) 


435.3 


5.6 


201.9 


1.1 


T Ba -Bb (MHJD) 


52941 


20 


53013.5 


7.8 


eBa-Bb 





(Fixed) 


0.67 


0.24 


acOL (/^as) 


221 


40 


217 


57 


*Ba-Bb (deg) 


87.2 


4.6 


84.3 


3.4 


WBa-Bb (deg) 





(Fixed) 


230 


23 


^Ba-Bb,l (deg) 


128.3 


4.0 


124.2 


3.9 


X 2 and dof 


1551.4 


1328 


1549.1 


1326 



Note. — Possible orbits for HD 196524, in the Campbell basis with 
lcr uncertainties. 

found that a 10 M® core could grow in ~ 10 Myr, but 
the core always formed at a distance of 1.5 AU, rather 
than at the observed 2.1 AU. Protoplanetary disks are 
seldom observed to survive for ~ 10 Myr around single 
young stars, much less binary stars, making core accre- 
tion appear to be an unlikely formation mechanism for 
gas giants in relatively close binary star systems. 

The alternati ve giant pl a net formation mechanism is 
disk instability ( |Boss|1997[ ). |Nelson| p000[ ) modeled disk 
instabilities in an equal-mass binary system with semi- 
major axis a = 50 AU and eccentricity e = 0.3, but found 
that the disks became too hot to fra gment into ga s giant 



protoplancts. On the other hand, Mayer et al. (20051 
found that disk instabilities could form gas giant planets 
in binary systems with e = 0.14 and a = 116 AU, but 
with a = 58 AU, whether fragmentation occurred or not 
depended on the protoplan etary disk m asses and the as- 
sumed disk cooling rates. Boss (2006) found that disk 
instabilities could lead to giant protoplanet formation in 
binary systems with semimaj or axes of 50 or 100 AU and 
eccentricities of 0.25 and 0.5. Mayer et al.| ( |2007| tried to 
reconcile these disparate results for disk instability, but 
could only conclude that given the problems with core 
accretion and the observational fact that gas giants ex- 
ist in binary star systems, disk instability remained as 
a possible formation mechanism for such planetary sys- 
tems. 

The candidate substellar companions discovered by 
PHASES require continued observations by other meth- 
ods for confirmation. Because PTI ceased operations in 
2008, acquiring new PHASES observations will not be 
possible. It is unlikely any existing northern hemisphere 
long baseline interferometers have the stable astrometric 
baselines required for differential astrometry, though the 
Navy Prototype Optical Interferometer may be a can- 
didate site. However, it would be better if an indepen- 



dent metho d could be used. Recent work by |Helminiak 



et all (120091) shows 40-1000 /ias precision astrometry us 



mg adaptive optics (AO) on large telescopes, while La- 
zorenko et al. (2007 2009) show similar precisions with- 



out AU on larger fields, which in principle might be ap- 
plied to AO images capable of resolving the binaries. 
The PHASES candidate systems should be high prior- 
ity targets for those observing programs; it is likely their 



precisions are sufficient to confirm or reject most of the 
candidate companions. The SIM-Lite Astrom etric Ob- 
servatory ( Shao et al.||1995 Unwin et al.||2008 ) will also 
be capable of confirming these companions and identify- 
ing additional systems. 
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ment. Without this outstanding engineering effort to 
produce a solid foundation, advanced phase-referencing 
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